The Observational Implications of Loop Quantum Cosmology 
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In this paper we consider realistic model of inflation embedded in the framework of loop quantum 
cosmology. Phase of inflation is preceded here by the phase of a quantum bounce. We show how 
parameters of inflation depend on the initial conditions established in the contracting, pre-bounce 
phase. Our investigations indicate that phase of the bounce easily sets proper initial conditions 
for the inflation. Subsequently we study observational effects that might arise due to the quantum 
gravitational modifications. We perform preliminary observational constraints for the Barbero- 
Immirzi parameter 7, critical density p c and parameter A. In the next step we study effects on power 
spectrum of perturbations. We calculate spectrum of perturbations from the bounce and from the 
joined bounce+inflation phase. Based on these studies we indicate possible way to relate quantum 
cosmological models with the astronomical observations. Using the Sachs- Wolfe approximation we 
calculate spectrum of the super-horizontal CMB anisotropics. We show that quantum cosmological 
effects can, in the natural way, explain suppression of the low CMB multipoles. We show that 
fine-tuning is not required here and model is consistent with observations. We also analyse other 
possible probes of the quantum cosmologies and discuss perspectives of their implementation. 



I. INTRODUCTION 



A main obstacle in formulating quantum theory of 
gravitational interactions is the lack of any empirical clue. 
Here the problem is that quantum gravity effects are ex- 
pected to be significant at the energies approaching 10 19 
GeV (the Planck scale). With the present generation of 
accelerators, energies up to 10 3 GeV can be reached what 
is sixteen orders of magnitude below the Planck scale. 
Therefore direct experimental investigation of quantum 
gravity effects becomes inaccessible. In other words it is 
like probing atomic structure with Earth size resolution 
devices. This suggest that alternative methods of inves- 
tigation ought to be taken into account. One, perhaps 
most promising, possibility of escape from this impasse 
are indirect methods. In this paper we consider one par- 
ticular type of indirect probing of quantum gravitational 
effects, which based on cosmological observations. 

In order to perform any quantitative predictions, 
mathematical model of the given process has to be 
known. The process considered here is behaviour of the 
universe in the Planck Epoch. In this epoch, evolution of 
the universe is determined by the quantum gravitational 
effects. In our considerations this phase is described 
by loop quantum cosmology (LQC) LQC base on 
nonperturbative approach to quantise gravity called loop 
quantum gravity (LQG) @]. Starting point in formulat- 
ing both LQG and LQC is paramctrisation of the phase 
space of the gravitational field by ST/ (2) connection and 
by its conjugated momenta. These canonical fields are 
so called Ashtekar variables (A = A l a r l dx a ,E = Efr l d a ) 
which take value in su(2) and su(2)* algebras respectively 



and fulfil the Poisson bracket 



{4(x),R 6 (y)}= 7K ^ 3 )(x-y) 



(1) 



where n = 8?rG and 7 is Barbero-Immirzi parameter. 
Parameter 7 is unknown constant of the theory. How- 
ever, because 7 is related with a black hole entropy, its 
value can be recovered from comparison with Hawking- 
Bekenstein formula Sbh = ttt- A. In our considerations 

4( pl 

we use value 7 = 7m = 0.2375 calculated by Meissner in 
Rcf. 0). 

Loop quantum cosmology can be considered on the 
two levels: the first is purely quantum approach and 
the second is semi-classical, effective framework. The 
first approach is more complete with respect to effective 
approach. However the semi-classical approach is more 
useful in relating quantum cosmological effects with clas- 
sical physics. Moreover, main features of the complete 
approach are sufficiently well reproduced by the effec- 
tive approach. Because our aim here is to relate effects 
of LQC with subsequent classical evolution, the semi- 
classical approach will be more adequate. Therefore in 
all the considerations performed in this paper we base on 
semi-classical LQC. 

In the cosmological applications, canonical variables 
(A, E) can be split for the homogeneous and perturbation 
parts: 



A 1 — A 1 

-fin 



SAl 



and E] = E] 



5E° 



(2) 
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In this paper, cosmological background is described by 
the flat FRW metric, then A l a = pS l a and £* = -fkS^. 
Here k and p are new canonical variables which fulfil the 
Poisson bracket {k,p} = The parameter Vo is the 
fiducial cell which regularise integration over the infinite 
spatial part. The p variable can be expressed in terms of 
the scale factor, p = a 2 and k = p/2p. Having canonical 
variables one can introduce Hamiltonian. The Hamilto- 
nian can be also decomposed for the homogeneous and 
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^phc„_ 



<5#g hen . Here the 



the background parts, Hq' 
upper index symbolise that classical Hamiltonian con 
tains additional phcnomcnological quantum correction. 
The lower index indicate that the gravitational part is 
considered. However in the realistic models, the matter 
Hamiltonian also contribute, then H phen = Hq 1011 + H m . 
In the considered models there is no quantum correc- 
tions to the matter Hamiltonian, what is indicated by the 
lack of the upper index. The matter Hamiltonian can be 
also decomposed for the homogeneous and perturbations 
parts. We will analyse such a case in this paper. 

In this paper we consider effective LQC with a scalar 
field. In particular we concentrate our attention on the 
model with a massive potential. In this case it will 
be possible to investigate realisation of the inflationary 
phase. In this approach, parameters of inflation arc 
dependent on the previous quantum cosmological evo- 
lution. We will consider perturbations in the emerg- 
ing bounce+inflation scenario. This will allow us to re- 
late quantum cosmological effects with classical pertur- 
bations. This finally let us to study LQC modifications 
of the CMB anisotropy. 

The organisation of the text is the following. In Sec. 
[II] we introduce concept of cosmic bounce in the frame- 
work of the effective LQC. Subsequently in Sec. IIIII wc 
construct a model of inflation in the applied framework. 
We set initial conditions in the contracting phase and 
study how parameters of inflation vary with them. We 
show that phase of bounce can easily set proper ini- 
tial conditions for the subsequent phase of inflation. In 
the next step, in Sec. HVi wc discuss perturbations in 
the considered model. We restrict ourselves to the fluc- 
tuations of the scalar field. We calculate spectrum of 
the perturbations from the bounce and from the joined 
bounce+inflation phase. These results can be applied 
in calculating spectrum of the CMB anisotropies. In 
Sec. [V] based on Sachs- Wolfe approximation, we calcu- 
late spectrum of temperature anisotropies in CMB. We 
show that phase of bounce can lead to suppression of 
the low multipoles in the spectrum of CMB anisotropics. 
Subsequently, in Sec. IVI1 wc discuss other kinds of the 
observational probes of LQC. Finally in Sec. IVIII we sum- 
marise our results. 



where A = 2y/3Trjlp l . In our investigations we use this 
particular expression for p. The choice of A is moti- 
vated by the existence of the gap in the spectrum of area 
operator in LQC However, A is the result of the kine- 
matic sector of LQG and extrapolation of this result to 
LQC is an assumption. This issue is discussed in Rcf. 
[J, The problem here is that the relation between 
LQC and LQG is not fully understood. In particular 
it should be possible to derive LQC directly from LQG 
and then problem of ambiguities in LQC should be over- 
came. Recently one promising step has been done to- 
wards to derive LQC from LQG. In their work, Rovelli 
and Vidotto show how LQC can be derived in the spin 
networks formalism. Another possibility is those pre- 
sented in where p, = X/y/p and A is some unknown 
constant unrelated with A. From this point of view A is 
some phcnomcnological parameter which should be de- 
termined from the observations rather from the theory. 

Taking Hamilton equation p = {p, H m + Hq 1011 } to- 
gether with the scalar constraint H m + Hq 1011 = one 
can derive modified Fricdmann equation 



where 



H 2 := 



Pc 



1 dp 
2p~dt 



P 1- 



V3 



16tt 2 7 3 ^ 



= 0.82 rri 



pi 



(4) 



(5) 



pi 



is critical energy density. As one can find from Eq. (j4|) 
the physical solutions are allowed only for p < p c . The 
p = p c is a turning point, commonly called a bounce. 
Moreover, while p — > the classical dynamics is recov- 
ered. One can also find that maximal value of the Hubble 
factor defined in Eq. (U]) is reached for p max = ^. The 
maximal value of the Hubble factor is -ff^ax = j^Pc- In 
Fig. [T] we show typical symmetric bounce obtained for 
the model with a free scalar field. 



II. BACKGROUND DYNAMICS 

In our considerations, Hamiltonian of the gravitational 
homogeneous part is given by 



H, 



phen 
G 



3Vb f= ( sin pjk 



-VP 



in 



(3) 



The crucial element of this Hamiltonian is the factor p. 
This parameter contains details of the quantum modifi- 
cations and when p — > 0, the classical limit is recovered. 
The main ambiguity in LQC comes from the choice of 

p. The mostly used form of p is this given by p = J 4 



FIG. I: Evolution of variable p and Hubble factor H in the 
bouncing universe (thick lines) with a free scalar field. Dashed 
lines represents classical evolution. Dots represents the points 

(t±,T-ffm«) = (±7K£.TV^j 
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III. QUANTUM BOUNCE AND INFLATION 



A. Analytical approximations 



In this section we are going to construct realistic model 
of inflation embedded in the framework of effective loop 
quantum cosmology. Qualitative studies of inflation in 
LQC have been already performed in Ref. 0]. Issue of 
inflation in LQC has been raised also in Ref. [{| Q . In 
our studies we model phase of inflation with the massive 
scalar field. Since we consider flat FRW model the equa- 
tion for the homogeneous component of the field p holds 
the classical form 



d 2 p 

dt 2 



3H ^ + dv 

at dp 







(6) 



where massive potential is given by 

v(p) = ^V. 

Energy density of the considered homogeneous scalar 
field is 



(7) 



(8) 



Dynamics of the model is governed by the set of equa- 
tions 



dH 
~dt 
dp 

~dt 

dip 
~dt 
dP 

dt 



P 2 



-£ = -p. 



— = -3HP- 



p2 

2 



V(<p) 



1 



dp 



(9) 

(10) 
(11) 
(12) 



Phase space of this system has been studied in Ref. 
Q. Analogous dynamics for the closed FRW model has 
been studied recently in Ref. [TFJ| . 

In our consideration we are going to restrict ourselves 
to the subset of initial conditions. We consider initial 
condition in the pre-bounce stage at the time to- We 
make very general assumption that at this arbitrary time, 
field is placed in the bottom of the potential, therefore 
ip(to) = 0. Since p < p c we obtain restriction for P at to, 
namely |P(to)| < y/2p c - Taking particular value of P(to) 
we can compute value of the Hubble factor 



H(t ) = -d- 




(13) 



In the dynamical system defined by Eq. [51 1121 one can 
distinguish subsystem (H, ip, P) whose evolution does 
not depend on p. In this subsystem, initial conditions 
are specified by the value of P(to), because (p(to) = and 
H(to) is given by Eq. [12J Initial value of p can be set 
arbitrary since only changes of p have physical meaning. 

In the subsequent part of this section we will show how 
parameters of inflation depend on the choice of P(to) and 

TO. 



Dynamics of the considered model is nonlinear and 
cannot be traced analytically. However we can distin- 
guish two regions where approximated analytical solu- 
tions can be found. The first is the phase of contraction 
and the second is the phase of slow-roll inflation. In 
this first phase, field oscillate in the bottom of the po- 
tential well. Therefore value of p = is reached many 
times, what motivate our choice of the initial condition 
p{to) = 0. The oscillations are amplified when the mo- 
ment of the bounce is approached. In this regime evolu- 
tion of the field is approximated by 



(pit) = C 



cos [m(t — 1 )} 



P 



3/4 



(14) 



In the subsequent phase of the slow-roll inflation, evolu- 
tion of the scalar field is approximated by 



pit) = p n 



-J. 



(15) 



In Fig. [2] we show exemplary evolution of the scalar field 
for the considered model. We show also how approxi- 
mated solutions fit to the solution obtained numerically. 
In the contracting phase, scalar field follows the approxi- 



mt 



FIG. 2: Evolution of the field ip. Dashed lines represents 
analytical approximations. Here m = f0~ 4 mpi and <^ ma x = 
2.9m P i. 

mated solution given by Eq. Q3] Therefore energy density 



behaves like p 



This is equivalent with the case 



of the universe filled by the dust matter. Corresponding 
evolution of the parameter p is in this case given by 



Pit) 



'3k 



Cmt + p^ 4 



1/3 



(16) 



In the subsequent phase of the slow-roll inflation, the 
approximated solution is given by 



p(t) = pi exp 



16ttG / _ to 
— 5 — m ip max t 



(17) 
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In Fig. 


[3] we show exemplary evolution of the canoni- 


P(t )/m 


1 


to- 1 


10 -2 


10~ 3 


10~ 4 


10~ 5 


10~ 6 


cal variable p for the considered model. We show also 


1 


0.5 


0.8 


1.1 


1.4 


1.8 


2.1 


2.2 


how approximated solutions fit to the solution obtained 


to- 1 


0.9 


1.1 


1.5 


1.8 


2.2 


2.4 


2.7 


numerically. 


10~ 2 


0.7 


1.6 


1.8 


2.2 


2.5 


2.8 


3.0 






1(T 3 




1.3 


2.2 


2.5 


2.9 


3.2 


3.4 


V 


P 


HT 4 






2.0 


3.0 


3.2 


3.6 


4.0 


10 2 ' 












2.7 


3.7 


3.9 


4.2 
















3.4 


4.4 


4.5 


10 2 ' 


/ s 

i \ 
/ \ 


lO" 7 












4.1 


5.0 



1(1" - 



mt 



FIG. 3: Evolution of the variable p. Dashed lines represents 



analytical approximations. 
2.9m P i. 



Here m = 10 mpi and (^max = 



FABLE I: Values of ip ma x for the different -P(io) and m (all 
parameters in Planck units). 



change considerably Therefor no fine-tuning is required 
to obtain proper phase of inflation. Moreover it is worth 
to stress that phase of bounce indeed leads to the proper 
inflationary scenario. In the classical considerations the 
high value of </? max has to be assumed, while in the LQC 
inspired model, this value can be obtained naturally. 



B. Conditions for inflation 

When the field (p reaches a point of maximal displace- 
ment (/? max then turns back and consequently P(<^ ma x) = 
0. At this point energy of the field is given only by the 
potential part. Because the total energy density is re- 
stricted p < p c , we obtain the following constraint 

|*W| < — = 1-3^. (18) 
to m 

The parameter (p max is important since gives good 
characterisation of the inflation. Moreover based on its 
value one can express e-folding number as follows 

- 9 

7V~27r^2f^. (19) 

Based on this expression and Eq. [18] we obtain another 
bound 

" m P1 

In the bounds given by Eq. [TBI and Eq. [20] we have as- 
sumed value of p c given by Eq. O However these bounds 
can be seen also in the different way. Namely, having 
parameters of inflation one can restrict the value of p c . 

In the considered setup the value of (/? max depends only 
on P(to) and m. It is worth to study how the value of 
<Anax is sensitive on them. The results of our investiga- 
tion are shown in Tab. [TJ 

In this table we collected values of ^ max obtained for 
the different values of initial parameters. The main con- 
clusion coming from this data is that despite the substan- 
tial change of the parameters, the value of (/? max does not 



C. Constraining p c , 7 and A 

Now let us assume that parameters of inflation are 
given by to = 10 _6 topi and </J max = 3.4topi, what gives 
~ 73. These values comes partially from the CMB 
observations and partially from the the requirement of 
solving the horizon problem. Based on these values and 
from Eq. [15] we obtain the following constraint 

p c > 6- l(r 12 TO P1 . (21) 

This constraint is not very useful because is very week. 
However the point is just to show possibility of con- 
straining and indicate the presently available cosmolog- 
ical bounds. Lets us also examine restriction of the 
Barbero-Immirzi parameter. Taking Eq. [5] together with 
the constraint Eq. [21] we obtain 

7 < 1222. (22) 

Now one can also constraint the value of parameter 
A from the formulation presented in Ref. [fj, where 

/2 = X/y/p. Here A is phenomenological scale of the loop 

2 

quantisation (polymerisation). Now p c = and 
assuming 7 = 7m = 0.2375 we obtain 

A < 7 ■ 10%i. (23) 



IV. PERTURBATIONS 

Cosmological perturbations are essential clement in 
searching for the quantum gravitational signatures. It 
is because the super-horizontal modes of perturbations 
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can carry information from the inflationary and even pre- 
inflationary epoch. Another important issue is that gen- 
eration of perturbations can have quantum gravitational 
origin. 

In LQC, as we already mentioned in Introduction, the 
perturbations are introduced according to Eq. [5] Per- 
turbations (5 A, SE) can be split for the: 

• Scalar modes (coupled with a scalar field) 

This type of perturbations with LQG corrections 
was studied in Ref. [n], EH, EH- However until 
not only inverse volume corrections have been in- 
troduced systematically. Consistent introduction 
of the holonomy corrections to the scalar modes is 
still awaiting. 

• Vector modes 

This kind of perturbations are of the secondary im- 
portance in cosmology. It is because they are de- 
caying modes and cannot affect the CMB substan- 
tially. Nevertheless this type of perturbations was 
studied in the context of LQC in Ref. 

• Tensor modes (gravitational waves) 

Tensor modes in LQC were studied in numerous 
papers. In contrast to the other two types of per- 
turbations, in this case also phenomenological im- 
plications were studied. The effect of inverse vol- 
ume corrections were studied in Ref. [H, EE E3 
while effects of the holonomy corrections were in- 
vestigated in Ref. [H, QJ, EI M, HI- In these 
papers, creation of gravitational waves was consid- 
ered either during the phase of a bounce or during 
the phase of inflation. The next natural step here 
is to consider creation of the gravitational waves 
at the joined bounce+inflation phase considered in 

sec. nni 

I this section we will consider simplified model of per- 
turbations. Namely we will consider perturbations of the 
matter field only. The gravitational field is set to be ho- 
mogeneous. This is only idealisation, however many re- 
sults from these studies can be extrapolated to the case 
of scalar and tensor perturbations. 



A. Scalar field perturbations 

Hamiltonian of the scalar field is given by 



H„ = d 



1 K 



Vo \2 ^/IdctB 

1 EfE^d a pd b p 



2 ^\drtE\ 



+ y/\dctE\V(<p) . (24) 



Similarly like in the case of gravitational field, scalar field 
can be decomposed for homogeneous and perturbation 



parts 

ip = (f + 5tp tt v = tt v + Stt v . (25) 
Here homogeneous parts are defined as follows 

<p{t) = ±r [ rf 3 x^(x,i), (26) 

Jvo 

M*) = W f d 3 x^(x,i). (27) 

Equations of motion for the background and perturba- 
tion parts are given by 



<P = {<P,H<p} =P 3/2 ' 7T l p, 



Sip = {5(p,H tp }=p- 3 / 2 6ir 

= 

'a 



v — {Sir v ,H v } = (VpV Sip 



P 



dip 2 



(28) 
(29) 
(30) 

(31) 



Combining equations (f2"B"|) and (f2"9")) we obtain equation 
©. Variable Sip can be decomposed for the Fourier 
modes 



d 3 k u(?7,k) ik 



(32) 



Based on this decomposition and on the equations (|30j) 
and (|31[) we obtain equation 



dr/ 2 

where k 2 = k • k and 



(77, k) + [k 2 + m 2 cS ]u{r/,k) = 0, 



2 _ _d 2 V(p) 1 d 2 VP 



m cff = V- 



dip 2 y/p dr/ 2 



(33) 



(34) 



In order to describe properties of the perturbations it is 
useful to introduce quantity called power spectrum which 
is defined as follows 



k 3 , 
Vu = 2^ |U 



(35) 



In the following two subsections we compute this quan- 
tity for two quantum cosmological models. The first 
will be the model of the symmetric bounce with the free 
scalar field. The second will be the model with the joined 
bounce+inflation phase. 



B. Symmetric bounce 

As the first case we consider scalar field perturbations 
at the symmetric bounce. In the considered case the 
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field is free, V = 0. We set the initial conditions to be 
the Minkowski vacuum 



-ikrj 



2k 



(36) 



This approximation works however only for the sub- 
horizontal modes. With these initial conditions we solve 
the equation (|33| numerically. Based on these compu- 
tations we obtain the power spectrum of the field u in 
the post-bounce phase. We show the results in Fig. [4] 
In this figure, black straight line represent analytical ap- 




FIG. 4: Numerical power spectrum of the field u (green 
points). Black line represent the analytical spectrum given 
by Eq. f42]with Uo = 3mp! and 770 = 0.1. Dashed lines repre- 
sents UV and IR limits given by Eq. 03] and Eq. 03] 



proximation of the spectrum. In order to derive this ap- 
proximation we assume 



2k 



2k 



(37) 



Here the relation \ak\ 2 — \f3k\ 2 = 1 holds, as a conse- 
quence of the normalisation condition. Now we base on 
the integral representation 

1 fi 

u(t), k) = win (77, k) + - / dr\ U(r] ) sin k(r] — 77 k) 
k J -00 

(38) 

of the Eq. 1331 where U(j]) = —m 2 s (r]). Solving this equa- 
tion in the first order of perturbative expansion we com- 
pute values of and f3u- Performing approximation 
U(rj) = U @{r) - t]o)Q(t]o - rj) where 



-m 2 oS (t = 0) = - (27f vPc ) 2/3 



(39) 



wc find 



/OO r r 

dr 1 U(r 1 ) = l + i^, (40) 



Pk ~ -2k 



drjU (rj)e 



— 2ikrj 



l ^Lsm(2k Vo ). (41) 



Now with use of definition of the power spectrum we 
obtain 



U (sin [2fc?7o] 2 U + 4k 2 U Vo + 4fc sin [ 2fcr ?o] (k sin[2r]k] - cos[2r7fc]r7 ??o)) 



16nk 2 



(42) 



In the UV and IR limits, the power spectrum given by 
Eq. behaves like 



k 

27T, 

2tt 



(43) 



1 + iUoVoV + 4(7 2 % 2 77 2 ) . (44) 



The term in the second bracket in Eq. [44] is the difference 
between the UV and IR slopes in Fig. |4"21 

The spectrum obtained in this subsection is similar 
to this obtained in the case of gravitational waves in 
Ref. [Isj . The difference is that in the case of the 
scalar field the effective mass m 2 s is negative in vicin- 
ity of the bounce while in case of the gravitational waves 
with holonomy corrections, the effective mass is positive 
function during the whole evolution. 



C. Bounce+Inflation model 

Phase of symmetric bounce is very idealised situation. 
More physically relevant is the joined bounce+inflation 
phase. In this subsection we show simple analytical 
model of perturbations in this phase. It will be a model 
of perturbations created in the scenario described in Sec. 
IIII1 In the contracting phase the sub-horizontal modes 
are given by Eq. 1361 The subsequent phase of inflation is 
approximated by de Sitter phase where evolution of the 
scale factor is given by a — — -grj- . In this phase modes of 
fluctuations are given by superposition of Bunch-Davies 
modes 



"out 



-ikrj 



1 



I 

krj 



1 



1 

kr] 



(45) 



where relation |afc| 2 — \(3k\ 2 = 1 holds. Performing match- 



ing conditions Ui n (?7i) = u out (r]i) and u' in (r)i) 
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we determinate cocfHcients 

1 - 2ikr}i - 2k 2 r/ 2 



2k 2 rj 2 



Pk = - 



— 2ikr)i 

2k 2 n 2 



(46) 

(47) 
(48) 



Based on this we can derive the power spectrum. We 
show plot of this spectrum in Fig. O 

<P Slf /(H Q /2n) 2 




FIG. 5: Power spectrum of the field Sip. Here rji — 
-1,-10, -100 (from right to left). 

The obtained spectrum is characterised by the suppres- 
sion for the low values of k. For the large k the spectrum 
holds the inflationary form. Another important feature 
are oscillations which are residue of the bouncing phase. 
We see that damping begins when — r]ik ~ 1. This cor- 
responds to the k on horizon scale at time rji. Based on 
this we define fe* = — ~ . At the time r)i a value of the 

scale factor is given by a, = — g * , therefore fc* = a* .Ho- 
Defining the length scale A* = a;/fc* we obtain A* = -g-. 
Today value of the A* is given by X^ao/ai where do is the 
present value of a scale factor. 

The similar power spectrum to this obtained here was 
derived also in Ref. [22, 01 ■ 



V. CMB ANISOTROPY 

As we have shown, there are two effects of the bounce 
phase on the primordial power spectrum: damping of 
the low energy modes and oscillations. This first effect is 
dominant and in this section we are going to investigate 
its impact of the CMB anisotropy. 



A. Sachs- Wolfe approximation 

Because we expect that effects of the bounce can affect 
super horizontal modes, the Sachs- Wolfe approximation 
can be used to study resulting spectrum of CMB. In this 



approximation sub-horizontal evolution of the primordial 
plasma is neglected since it does not affect the considered 
modes. Expression for the CMB multipoles is given by 



C, 



4?r 
25 



dk 



V n (k)jf(kD, 



(49) 



where D+ = 1.4 • 10 4 Mpc is distance to last scattering 
shell. Moreover Vn(k) = ^|^-k| 2 where TZ = -f and v 
is Mukhanov variable which fulfils equation 



k z - — 



and where z 



H ' 



We see that while approximation z / ' z ~ a / 'a holds, 
then evolution of v and u variables are the same. The va- 
lidity of this approximation was indicated in Ref. 0, [23| . 
The assumption we made here is lack of the quantum 
modification to the equation for the v variable. However 
we do not have a right equation for the scalar modes in 
presence of holonomy corrections. In case of the tensor 
modes, it was shown in Ref. [IH that corrections to the 
mode equation do not change the spectrum significantly. 
It was shown that the shape of the spectrum is deter- 
mined mainly by the background evolution. Therefore 
we assume here that the spectrum of the scalar pertur- 
bations is not affected significantly by the holonomy cor- 
rections to the mode equation. Then spectrum from the 
joined bounce and inflation phase should have generic 
form derived in Ref. IIVI In order to build the analytic 
model we can average the spectrum over the substandard 
oscillations. Then the power spectrum from the joined 
bouncc+inflation is approximated by 



V n {k) = A^Q{k - h) + A%e{h - k) 



(50) 



Based on this spectrum one can derive analytical for- 
mula for the spectrum of the low multipoles of the CMB 
anisotropy. Instead of the variable C\ it is convenient to 
consider variable 



(51) 



It is motivated by the fact that for the scale invariant 
Harrison-Zeldovich power spectrum, this quantity holds 
constant value. Performing integral (|49p with the spec- 
trum (1501) we obtain 



where 



and 



^inflation _j_ ^bounce 



A 2 

^inflation 



25 



(52) 



(53) 
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.bounce = [ iai , 3(1 + + + 3/2; 2 + i, 21 + 2; - (1 + l)iF 2 (l; I + 3/2, 2? + 2; -aft] 



Here we have introduced parameter x„ — k^D*. In Fig. [5] 
we show Ci spectrum with parameter A\ — 2.6 ■ 10 -9 set 
to fit the CMB data and T CM b = 2.726 K. We see that 




FIG. 6: Spectrum of CMB anisotropy. 

effects of the bounce lead to suppression of the low mul- 
tipoles in CMB, what is favoured observationally. This 
possibility was indicated earlier in Rcf . 0, [23[ . 

Value of parameter A\ can be calculated from the 
slow-roll inflation model 



A 2 



1 



-V 

271"/ 



where e is slow-roll parameter 



1 

2k 



1 dV 

Vd£ 



'pi 



4ir (p 2 



(55) 



(56) 



Based on this we can calculate value of the Hubble factor 
at the beginning of inflation 



H = 



'PI 



= 3.8 • lCr 5 m P i 



(57) 



where in the last equality we assumed <p max = 3.4mpi. 

Important limitation of the method based on the low 
multipolcs in CMB comes from the so called cosmic vari- 
ance. It is purely statistical effect which is significant at 
the horizontal scales. In case of CMB this corresponds to 
the low multipoles regime. Relative uncertainly coming 
from the cosmic variance is given by 



AC; 



21 + 1 



(58) 



Taking for example I = 2 we obtain relative uncertainty 
equal 0.63. Therefore outcome of measurement is com- 
parable with its uncertainty. This effect cannot be re- 
moved and impose substantial limitation on our approach 
to constraint quantum cosmological models. 



B. Discussion 

In the previous subsection we showed that bounce can 
leads to observed suppression of the low CMB multipoles. 
Now we would like to discuss whether this scenario is re- 
alistic and does non require fine-tuning. At the beginning 
we would like to however mention one important adjust- 
ment. Namely observed scale of cut-off in CMB spectrum 
overlaps with present size of cosmic horizon. This prop- 
erty was indicated in Rcf. '22^. Therefore there is an 
intriguing possibility that observed cut-off is due to un- 
known evolution of the super-horizontal modes. There- 
fore it is not related with the pre- inflationary dynamics. 
Such an explanation seems to be more likely. However 
model of this super-horizontal damping does not exist 
yet. 

It is indicated by the observations that present scale 
of cut-off A* (t ) is comparable with A* (to) ~ D*. 
Therefore 



a* 



(1 + Zdcc) 



(59) 



Because A*(tj) = 1/H , taking Eq. IBT1 and Eq. IT91 we 
obtain 



2Ne 



2N 



where 



2Dl{2vYml l Al ( T d 



(1 + Zdec) 2 



GUT 



(60) 



(61) 



It is worth to note that Eq. [60] has a form of the Lam- 
bert equation W{z)e w ^ = z which define Lambert W 
function, therefore N = §W(£). In order to determi- 
nate parameter £ we take za cc ~ 1070, Td ec — 0.2 cV, 
T GUT ~ 10 14 GcV and A\ = 2.6 • 10" 9 . Based on this 
we obtain £ = 5.1 • 10 62 and subsequently N = 69.7. 
Now one can determinate second independent parameter 
of inflation e.g. m. We can easily derive equation 



3 ,o 27r 

— AZ, — mpi. 

2 n N 



(62) 



which gives m = 5.6T0 _6 mpi. As we see the obtained pa- 
rameters of inflation are fully consistent with these usu- 
ally considered. The model is in the full agreement with 
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the present observational facts. Moreover any fine-tuning 
is not needed to explain suppression of the low multipolcs 
by the quantum cosmological effects. 



VI. OTHER OBSERVATIONAL PROBES OF 
LQC 

Beside the effect of suppression of the low multipolcs 
also other potential probes of quantum cosmologies are 
available. In this section we review four possible ap- 
proaches. 



A. Polarisation of CMB 

In Sec. [V] we have shown how the quantum cosmo- 
logical effects can be related with the spectrum of CMB 
anisotropy. This method gives us one possible approach 
to constraint quantum cosmologies. However observa- 
tions of CMB brings us much more information than only 
anisotropies of temperature. Another important mea- 
sured quantity is polarisation of CMB radiation. This 
polarisation can be described by the spectrum, which de- 
pends on the primordial perturbations. While E-typc po- 
larisation depends on the both scalar and tensor compo- 
nents of perturbations, the B-type polarisation depends 
only on the tensor component. Spectrum of the E-type is 
already observed and can be used e.g. to constraint cut- 
off of the power spectrum. Recently such a study were 
performed in Ref. [23 ]. Based on the joined anisotropy 
and polarisation data it was shown that scale of cut-off 
in the power spectrum is C = ir ^ c _ t < 4.2 at 95% 
confidence level. While the constraint based only on the 
polarisation gives C < 5.2. This result show that po- 
larisation is a good tool to constraint the cut-off in the 
power spectrum. In particular while the considerations 
based only on the CMB anisotropies indicate a cut-off, 
the joined anisotropy and polarisation data indicate limit 
on the cut-off. This is crucial while constraining quantum 
cosmologies based on cut-offs resulting them. 

The second, B-type polarisation, still remains unreach- 
able observationally. However there is presently a huge 
effort to detect it. In particular mission like PLANCK 
aims to detect this type of polarisation. If the B-type po- 
larisation will be measured then amplitude of the tensor 
power spectrum can be determined. In case of slow-roll 
inflation this amplitude is given by 

Because H 2 — §p, the measurements of the B-type po- 
larisation enable us to determine energy scale of inflation. 
Therefore absolute values of parameters -/V and m can be 
found. This give us also restriction on the pre-bouncc 
initial condition, in particular on P(to). 



B. Non-Gaussianity 

When different modes of perturbations interact one 
another then field becomes non-Gaussian. This interac- 
tion can be produced by the potential of the scalar field, 
higher order corrections in the perturbative expansion or 
by the quantum gravitational effects. 

In case of the Gaussian field, all its statistical prop- 
erties are fully described by the two point function 

(PkiVk,) = ^ (3) (ki - ^)j£-V v (k) where T v (k) is power 
spectrum given by Eq. 1351 In order to describe non- 
Gaussian effects it is necessary to consider higher or- 
der correlation functions. The first contribution of non- 
linearity comes in tree-point function 

(VkiPkaVk.) =<5 (3 )(ki+k 2 +k 3 )P3(ki,k 2 ,k 3 ) (64) 

where _P3(ki, k 2 , ka) is called bispectrum. In case of the 
Gaussian field, this spectrum is equal zero. 

The primordial non-Gaussianity, if present, could af- 
fect also the spectrum of CMB anisotropies leading to 
its non-Gaussianity. Present observations indicate how- 
ever that CMB spectrum is nearly Gaussian. This give 
us constraint on the cosmological models with a huge 
amount of nonlinearity. In particular, based on these 
observations, some quantum cosmological models can be 
constrained or even excluded. For example preliminary 
studies on non-Gaussianity in LQC were performed in 
Ref. J25|. In these studies, non-Gaussianity is produced 
by the specific scalar field potential not by the quantum 
gravity effects itself. However, this model gives an exam- 
ple of non-Gaussianity production in the bouncing uni- 
verse. It was shown that in this model, non-Gaussianity 
is produced in the vicinity of the points t + and t_, where 
Hubble factor reach its maximal value. Another example 
of non-Gaussianity production in the bouncing cosmol- 
ogy can be found in Ref. [26|. In this paper production 
of the non-Gaussianity at the matter bounce is consid- 
ered and indicate that this form of non-Gaussianity can 
be potentially distinguished from this produced during 
the inflationary phase. 



C. Transplanckian modes 

When the length of the mode of perturbation ap- 
proaching the Planck scale then semi-classical approxi- 
mation fails. The notion of continuous wave is missed 
and and fully quantum gravitational considerations have 
to be applied. Therefore these so called transplanckian 
modes cannot be studied with use of Quantum Field The- 
ory on Curved Spaces as it was done in the present paper. 
More adequate would be application of Quantum Field 
Theory on Quantum Spaces as this studied in Ref. [27| . 

In Fig. [7] we show evolution of the different length 
scales in the bouncing universe. We see that modes with 
A > l c can be described by the semi-classical approxi- 
mation. It is exactly the case considered in this paper. 
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FIG. 7: Evolution of the different length scales in the bounc- 
ing universe. 

However when A < l c then new formulation has to be ap- 
plied. This would lead to potentially new effects. How- 
ever these transplanckian modes can decay before cross- 
ing the horizon during the inflation. Then any signature 
of the quantum gravity effects can be lost. However in- 
vestigations as this performed in Ref. [H| suggest that 
effects of transplanckian modes can lead to potentially 
observational effects. 



D. Large scale structures 

If the quantum cosmologies can give imprints on CMB, 
then some of these effects could be seen also in the sub- 
sequent structures. The region of the low multipoles cor- 
responds now to the largest visible distances in the Uni- 
verse. Gravitational structures on these scales are called 
large scale structure (LSS). Therefore observations of the 
large scale structures are complementary to the observa- 
tions of the low multipoles in CMB. These both methods 
were already applied to investigate effects of the bounc- 
ing cosmological scenario in Ref. [29j . In this paper au- 
thors predict oscillations in the power spectrum on the 
horizontal scales due to the bounce. However available 
observational data from e.g. SDSS arc still not sufficient 
to probe these effects. 

VII. SUMMARY 

In this paper we have shown possible way to relate 
quantum cosmological models with astronomical obser- 
vations. This method can be used to constraint models of 
the universe in the Planck Epoch. In particular we con- 
strained Barmcro-Immirzi parameter 7, critical density 
p c and parameter A. However, present constraints are 
still very week and ambiguous. The part of this ambigu- 
ity can be removed by fixing the parameters of inflation 



(N and m). Parameter N is today fixed by the require- 
ment of solving the horizon problem. Then value of m 
is determined from the amplitude of CMB anisotropies. 
However problem of horizon do not appear in the bounc- 
ing cosmologies. Therefore at present we have fixed only 
relation between A an m but not the absolute values. 
In order to fix one of these parameters, other observable 
has to be measured. The most promising is amplitude 
of tensor perturbations. These perturbations produce B- 
type polarisation in CMB. Therefore, if this polarisation 
would be measured then parameters N and m can be 
fixed. This will enable us to perform less ambiguous con- 
straint on the LQC parameters, in particular on p c . 

We have shown that phase of a bounce set proper 
initial conditions for inflation. Moreover parameters 
of inflation depends only logarithmically on the pre- 
bounce initial conditions. Subsequently we considered 
model of perturbations at the bounce and at the joined 
bounce+inflation phase. We showed that this second 
model can give explanation of the suppression of the low 
multipoles in CMB. Moreover we have indicated that 
model is fully consistent and scale of the cut-off agree 
with the present size of horizon. This indicate that fine- 
tuning is not required to produce suppression on the hor- 
izontal scales. 

Beside the cosmological approach to constraint quan- 
tum gravity effects also other indirect methods are avail- 
able. In particular astrophysical measurements of the 
Lorentz symmetry violation. It is in principle possible 
to derive quantitative predictions about this effect di- 
rectly from LQG. However process of derivation require 
construction of the semi-classical states, where unknown 
phenomenological parameters appear. Since their val- 
ues are unknown, predictive power of this approach is 
marginal. If these difficulties would be removed then this 
method can be complementary with the cosmological ap- 
proach. Also the quantum effects on the gravitational 
collapses gives possible way to constraint the LQG. Here 
the results are less ambiguous, however it is harder to 
relate them with some astrophysical data. 

As we see the difficulties lies here on both theoreti- 
cal and empirical side. Without knowing the relation 
between the LQC and LQG we can treat parameters of 
this first rather as phenomenological. If these difficulties 
would not be overcome it will be hard to perform com- 
plementary constraints of the same parameters, based on 
different methods, e.g. observations of CMB and gamma 
ray bursts. 
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